Superconductivity in strongly repulsive fermions: the role of kinetic-energy frustration 
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We discuss a physical mechanism of a non-BCS nature which can stabihze a superconducting state 
in a strongly repulsive electronic system. By considering the two-dimensional Hubbard model with 
spatially modulated electron hoppings, we demonstrate how kinetic-energy frustration can lead to 
robust d-wave superconductivity at arbitrarily large on-site repulsion. This phenomenon should be 
observable in experiments using fermionic atoms, e.g. "^''K, in specially prepared optical lattices. 
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Introduction.- One of the long-standing fundamental 
questions in condensed-matter physics is whether it is 
possible to realize a superconducting (SC) state in a 
system consisting only of electrons subject to a strong 
Coulomb repulsion, and if so, what is the minimal set of 
necessary physical assumptions. An early attempt to pro- 
vide an answer was made by Kohn and Luttinger flj , who 
proposed a weak-coupling BCS-like mechanism. While 
' their idea was never confirmed experimentally, there ex- 
ist numerous strongly correlated systems whose SC be- 
havior occurs without any obvious pairing glue, such as 
phonons, between the electrons. Examples are high-T^ 
cuprates and heavy fermion compounds. The current 
consensus is that superconductivity in these materials 
has an unconventional, i.e. non-BCS, character Q. Un- 
derstanding the microscopic origin of this intriguing phe- 
nomenon remains a challenge. Here we address the above 
question by performing a controlled derivation of the SC 
ground state (GS) for a strongly-repulsive Hubbard model 
with spatially modulated transfer integrals. 
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FIG. 1: The tetrahedral lattice topology. The interplaque- 
tte hopping amplitudes are: NN ti (solid lines) and NNN T2 
(dashed lines). We only consider the case T2 ^ n- 



One possible way of stabilizing a Cooper pair conden- 
sate in a repulsive system is to introduce microscopic 
inhomogeneities. Indeed, the nanoscale spin and charge 
modulations, observed in scattering [3], ARPES (3] and 
STM [Ij experiments, seem to be ubiquitous in high- Tc 
materials [6( and often accompany the emergence of the 
SC state. Theoretically it has been argued that these in- 
homogeneities are quite relevant for the superconductiv- 
ity and, in fact, seem to assist the Cooper pairing. 
This was demonstrated in 0, 0| by using exact diagonal- 
ization of strongly interacting models in finite lattices. 
In Ref. [l3 the authors studied the Hubbard model on a 
checkerboard lattice, composed of weakly coupled 2x2 
plaquettes, and showed that the SC phase can be stabi- 
lized in a relatively narrow interval of the on-site repul- 
sion U. Earlier, a similar problem was considered in 11]. 
Another ingredient, whose importance for superconduc- 
tivity was largely overlooked, is the range of the transfer 
integrals beyond nearest- neighbors (NN). The next-NN 
(NNN) hopping, t' , was shown to enhance da;2_j,2-like 
pairing correlations in the t-t'-J model on finite clusters 
|12| . Physically, its main qualitative effect is the possible 



frustration of the kinetic-energy term: the smallest closed 
paths in the lattice are triangles instead of squares. 

In the present Letter we explicitly demonstrate how 
local kinetic-energy frustration can stabilize the SC state 
in a strongly repulsive two-dimensional Hubbard model. 
The lattice, on which the model is defined, is presented 
in Fig. [T] It consists of weakly-coupled tetrahedra, i.e. 
plaquettes with frustrated hoppings along the diagonals. 
We show that a c?a;2_j,2-wave SC phase exists for arbi- 
trarily strong repulsion U. In fact, the problem can be 
treated analytically in the strong-coupling regime. 

Our motivation to study this system is not purely aca- 
demic. Advances in experimental methods of preparation 
and manipulation of ultracold fermion atoms in optical 
lattices provide a controlled way of testing the above- 
mentioned theoretical ideas. For example, in recent ex- 
periments 14 1 the observation of a Mott state with 
'^'~'K atoms was reported. Moreover, an experiment aimed 
to find c?-wave superconductivity in a checkerboard Hub- 
bard model was proposed in Ref. [15I . 
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FIG. 2: Hole binding energy A(C/). The asymptotic behavior 
is: A{U t) ^ 2t^/U and A((7 < « /?,2t. Upper 
inset: critical value Uc{t') [A{Uc) = 0]. Lower inset: group 
Civ Numbers indicate plaquette vertices. The black square 
denotes four- fold axis C4, horizontal and vertical lines - pri- 
mary symmetry planes cr„, diagonals - secondary planes a'^. 



Model. - Let us consider the repulsive Hubbard model: 



H = - 



E 



(1) 



defined on the lattice, Fig. [1] in terms of fermionic (cre- 
ation) operators c|^. Here (ij) denotes links connect- 
ing sites i and j, a = {tii} is the electron spin, and 
"icr ~ Aa^iry- The amplitudes tij take four possible val- 
ues: (i) t for links (12), (13), (24) and (34); (h) t' for the 
diagonals (14) and (23); (iii) t\ for NN links, connecting 
two plaquettes; (iv) for NNN interplaquette links. 

We will consider the case ti_2 ^ t' , U, which allows 
for a controlled perturbative expansion of the Hamilto- 
nian ([1]). To demonstrate the existence of a robust SC 
phase, we derive a low-energy effective model, accurate 
to second order in T1.2. In general, this is doable only 
numerically. However, in the limit t, t' <C U, we can 
keep only lowest-order terms in tij/U, and thus provide 
a closed form for the effective Hamiltonian (EH). The 
stability of the Cooper pair condensate can be tuned by 
changing the ratio t'/t. There is an "optimal" value of 
this ratio, which ensures a finite energy gap (hole binding 
energy) between the plaquette states with one and two 
holes, for all finite U. 

Single-plaquette states. - The Hubbard Hamiltonian on 



a single plaquette can be diagonalized exactly [ll| by 
using representations of the crystallographic group C41, 
(see the lower inset in Fig. [2]). As a result, we can deter- 
mine the hole binding energy A = 2eo(3) — eo(2) — eo(4), 
where eo(-/Ve) is the GS for a given number of electrons 



iVg. Positive values of A correspond to binding of two 
holes. In general, A is positive only in a finite range of 
U. At some critical value Uc{t'), shown in the upper in- 
set of Fig. [21 it changes sign and remains negative as 
[/ — > 00. There is a special ratio, t'/t = 1, at which Uc 
diverges and A stays positive for any value of U (see the 
main panel of Fig. [2]). This results from the maximal 
frustration of the single-hole kinetic energy. The GS en- 
ergy for 4 electrons (zero holes), eo(4) — > for [/ — > 00 
because the particles cannot move. On the other hand, 
in this limit eo(2) = 2eo(3), which means that there is 
no kinetic-energy gain for creating two holes on different 
plaquettes; i.e., the single- hole kinetic energy is optimally 
frustrated. The exchange interaction J = At'^/U, that 
appears for finite t/U ^ 1, leads to pairing (A > 0) be- 
cause the magnetic configuration of two plaquettes with 
one hole in each of them is more frustrated than the con- 
figuration with two holes in the same plaquette. This 
leads to a positive value of A — J/2. 

From now on we will only consider the maximally frus- 
trated point t'/t = 1. Then, the symmetry group Q of 
the single-plaquette Hamiltonian is larger than C^y (sym- 
metry group for arbitrary t'/t), and contains all the in- 
dependent permutations of any pair of vertices of the 
plaquette. This symmetry translates into a GS degen- 
eracy at half-filling. There are two S'J7(2)-singlet states: 
one transforming as the identity representation of Cjy, 
Ai (s-wave), and the other - as Bi {d^-^ -y-i-'wsNe) [l6|. 
These states are connected by symmetry operations from 
the factor group Q/C^v The two-electron GS is also a 
singlet and belongs to the identity representation of Q. 
The Ne — i GS has S" = 1/2 and is six- fold degenerate. 

General expressions for these eigenstates are quite 
cumbersome. However, to the lowest order in t/U, 
we can consider only states without doubly occupied 
sites. Hence, we have the GS for = 2: \Vt2) = 
(l/2-\/3) summation extended over 

all links of a plaquette; and for Nf, = 4: = 
■^s.,d{s\'^s\^ ± s|2S34)|0). In these expressions s\^ is a 
singlet creation operator, sj^- = cj^^c^j^ 

empty state and Ms = —1/2, Nd = 
we introduce operators Pij, which permute sites i and 
j. In the basis {|i^|), l^^f)}, P12 and P13 have the form: 
-Pi2,i3 = -crV2±\/3cr^/2 with cr" [a = X, z) Pauh matri- 
ces. We will use this expression to determine symmetries 
of the effective model. 

Effective low-energy model - The low-energy spectrum 
of decoupled plaquettes has a gap A to single-hole {Ne — 
3 on each plaquette) states. Here we consider the ef- 
fect of finite hopping amplitudes ri^2 by assuming that 
^ ri,2 < A ~ t'^/U <^ t <^U. The second inequal- 
ity allows us to treat interplaquette hoppings perturba- 
tively. The fourth one allows us to exclude states with 
doubly occupied sites, i.e. use as a basis the states \Vl2) 
and jilj ). Finally, the third inequality constrains the 



c|'^c]^, |0) is the 
1/2 V^. Finally, 
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choice of the virtual states: only states that belong to 
the Nf, = 3 GS sextet contribute to lowest order. We will 
also assume that T2 ^ ti. 

The second-order EH can be symbolically written as: 



1 



En - i/(0) 



where H^^^ describes a set of noninteracting plaquettes 
in (dl, Eq is its GS energy, Hr denotes plaquette inter- 
actions, and V'^^^ is a projector onto the subspace with 
A'^e = 2 or 4 on each plaquette. Next, we associate the 
product of the two-electron plaquette GS with the vac- 
uum: |vac) = Yix \^2)x and each member of the four- 
electron GS doublet on plaquette x - with a hard-core 
boson: \fl2)x = ^'Icd^s)^, where a = s or d represents 
the pseudospin index. The algebra generated by b^a was 
discussed in Ref. |l7|. Thus, the effective low-energy the- 
ory, given by H^s, describes a system of two- flavor hard- 
core bosons, propagating in the coarse-grained plaquette 
lattice of Fig. [TJ In terms of these boson operators we 
have: 

HeS^ J2 t%{bLbvP + blf,bxo.)-tiJ2'^-°^+ (2) 

{xy),al3 x.a 

+ E yfp[i^-^^HJy0 + a-riy)bLbxp], 

{xy),aP 

where {xy) denotes NN plaquettes, are correspond- 
ing hopping amplitudes, - density-density and local 

spin-flip interactions, Uxa = bl^bxa and Ux = Uxs + rixd^ 
and /i is the chemical potential. Direct interactions be- 
tween pseudospins, like Heisenberg terms, are not present 
to lowest order in A/i and t/U. 

Some general properties of f^^ and can be estab- 
lished by symmetry arguments. First, the Hamiltonian 
([1]) is invariant under reflections in the planes which in- 
clude Ti links, e.g. the plane connecting sites 2 and 3 in 
Fig. [T] The states \Q.2) and are symmetric under 
this operation, while the d-wave state l^f) is antisym- 
metric. Consequently, the off-diagonal matrix elements 



of C and y/ff vanish: i^* 



'a/3- 



This result is independent of the assumptions made re- 
garding the relative magnitude of ri,2, U and t. 

Another observation concerns the diagonal elements of 
and V^^ in the special cases T2 = n and r2 = 0. In 
the first case, we consider the two plaquettes with num- 
bered sites, shown in Fig. [U and perform simultaneous 
permutations of vertices 1 o 2 on the left plaquette and 
1 O 3 on the right one. Each operation is a symmetry 
of the singic-plaquette Hamiltonian. Their combination 
amounts to interchanging the ti and T2 links, which is 
now a symmetry of the connecting Hamiltonian. Using 
the relation P12 131^^2) = \^2)^ it is easy to show that for 
Ti = T2: tf^ = -tf, and = V,f. In the case T2 = 0, 
when the plaquettes are connected by only one ri link. 



the second-order virtual hopping of an electron can only 
proceed through an intermediate state, whose energy is 
of order U . Therefore, in the approximation formulated 
above, f^^ must vanish. On the contrary, is not asso- 
ciated with the net electron transfer and remains finite. 

In general, a direct calculation yields the precise form 
of the coefficients t"^ and V°^: 



t' 



off 



„^ -- (r2/6A)diag{r,(2r,-|-l), -3r,}; (3) 
^ (rf /48A)diag{9 + 8r^ + IGr^, 9 + 24rl} 



T/eff 



with Vr = T2/T1. Clearly, in the two special cases, dis- 
cussed above - rr — 1 and ~ the EH ^ becomes pseu- 
dospin symmetric. The second case is irrelevant for the 
purposes of studying the SC state, while the first one, 
Tt- = 1, is quite instructive. Indeed, in this case we can 
use the Perron-Frobenius theorem to prove that there 
exists a pseudospin-polarized GS [11]. The Hamiltonian 
can then be written only in terms of spinless bosons, say 
bxd, and maps onto the spin- 1/2 XXZ model in a mag- 
netic field n, via the Matsubara-Matsuda transformation 
[1^ . The phase diagram of this model contains Neel, 
canted XY-antiferromagnetic and fully polarized states 
that are immediately identified with the density-wave 
(DW), Bose-Einstein condensate (BEC) of Cooper pairs, 
and Mott phases, respectively. The DW and BEC states 
are separated by a 1st order quantum phase transition. 

We do not expect the physics to change qualitatively 
for < r,- < 1. It is known that the usual mean-field 
approximation yields satisfactory results for Tt — I when 
compared to Monte-Carlo simulations [20|. Thus, we an- 
ticipate that the rest of the phase diagram, along the 
Vr axis, can be described within a simple variational ap- 
proach. We employ the method of 2l|, which includes 
short-range quantum fluctuations and, as a limit, con- 
tains the semiclassical spin-wave ansatz. The resulting 
phase diagram, obtained using 2x2 site clusters (in the 
plaquette lattice), is presented in Fig. |31 For any finite 
< r,- < 1 the system exhibits the same three phases, as 
in the case r,- = 1. The SC and DW phases are again sep- 
arated by a 1st order transition. The transition between 
SC and Mott phases is 2nd order. In the Mott state there 
is exactly one boson per site; i.e., the electron filling is 
1/2. In this phase the pseudospin polarization is unde- 
fined, as the Hamiltonian ^ becomes spin- independent. 
Interestingly, the DW phase is of an s-wave nature, due 
to the fact that the expectation value of the kinetic en- 
ergy vanishes, while density-density interactions favor the 
s-wave pseudospin polarization. 

The SC state has a dx2-y2-waYe symmetry. The struc- 
ture of the "Cooper pair" can be determined by observ- 
ing that bd = V/3 - (l/4)(si4 + S23)(II I) (I II " ^Oi), 
where |I I) = slgsLlO), IIH) = sl244|0) and V = S13 + 
S24 — S12 — S34 (see the lower inset of Fig. [2]). Hence, 
despite the apparent complexity of the SC phase, it can 
still be characterized by a familiar d-wave order parame- 
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ter Ad = (T)), shown in the inset of Fig. |3]for = 1. As 
Tr decreases, the height of the SC dome gradually dimin- 
ishes and disappears at =0. Thus, for any < r,- < 1 
there is an interval of /i where the SC phase is stabilized. 
This conclusion becomes rigorous in the dilute limit of 
particles or holes, by virtue of the inequality > |, 
vahd for r^- < 1 [see Eq. jS])]. 

Discussion. - Our phase diagram. Fig. [3l was obtained 
in the strong-coupling limit U ^ t, where one can derive 
the effective model of Eqs. ([2]), (|3]). The EH becomes 
increasingly complicated for intermediate couplings U ^ 
t, because of the large number of virtual transitions. In 
this regime, the existence of d-wave superconductivity 
in the nonfrustrated Hubbard model was argued in [lO| 
based on a first-order EH, treated within a mean-field 
approximation, and in the weak-coupling regime U <^ t 
in [22,]. Therefore, we expect the SC phase to persist for 
U ^ t in our frustrated case as well. However, regardless 
of the magnitude of U, the SC state is quite sensitive to 
the presence of longer-range repulsions. For instance, an 
interaction of the form J2 with Vij — V for all 

links within the plaquette, will suppress the local hole 
binding ifV ^Vc = O.lUt. For V < Vc the SC phase is 
stable only in a finite interval of U around U ^ It. 

Our theory highlights the importance of the kinetic- 
energy frustration for stabilizing the SC state. Locally, 
pairing competes against the kinetic energy and can be 
increased by frustrating the latter. This principle guides 
the choice of the elementary unit, e.g., tetrahedron. The 
connectivity of the lattice, built from these blocks is an- 
other essential ingredient. Here we used the lattice of 
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Fig. [T] to demonstrate the existence of the SC state in 
a physically transparent way. However, we also consid- 
ered the usual checkerboard lattice In this case the 
relation between coefficients in the EH is such that the 
phase-separated state can suppress superconductivity in 
a certain region of the phase diagram. The importance 
of the lattice topology is further illustrated by the case 
T2 = 0. Without the interplaquette hopping (jH]), the 
global phase coherence can be established only in higher 
orders in 1/C/, leading to a quite fragile SC state. 

Finally, we believe that the lattice of Fig. [T] can be 
realized using ideas of Refs. [iSl, |23|. Indeed, our effective 
strong-coupling model can be easily extended to the cur- 
rently experimentally realizable regime t' /t < 0.5 under 
the condition t <^U < Uc, which can still be fulfilled for 
t'/t = 0.5 because Uc ~ lit (see inset of Fig. [2]). The 
resulting phase diagram is qualitatively the same as the 
one shown in Fig. |31 Thus, results of the present Letter 
can be tested in future cold atom experiments. 
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FIG. 3: Low-energy phase diagram of the Hubbard model Eq. 
([l]). Phases are: s-wave density wave (DW) with wavevector 
(tt, tt); d-wave SC, which corresponds to a BEC of bd', Mott 
insulating (MI) phase with (ux) — 1. The DW - SC phase 
transition is 1st order; the SC - MI transition is 2nd order. 
The inset shows Ad ~ (bxd) for — 1. 
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